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Abstract. We study the fluctuations of the matrix entries of regular functions 
of Wigner random matrices in the Umit when the matrix size goes to infinity. 
In the case of the Gaussian ensembles (GOE and GUE) this problem was 
considered by A.Lytova and L.Pastur in 1281 . Our results are valid provided 
the off-diagonal matrix entries have finite fourth moment, the diagonal matrix 
entries have finite second moment, and the test functions have four continuous 
derivatives in a neighborhood of the support of the Wigner semicircle law. 
Moreover, if the marginal distributions satisfy the Poincare inequality our 
results are valid for Lipschitz continuous test functions. 



1. Introduction 

Let Xm — --j=Wn be a random Wigner real symmetric (Hermitian) matrix. In 
the real symmetric case, we assume that the ofF-diagonal entries 

(W^jv),fc, l<j<k<N, (1.1) 

are i.i.d. random variables with probability distribution /i, such that 

nWN)jk - 0, N{WN)jk = nWN)% - m4 < (X), l<]<k<N, (1.2) 

where denotes the mathematical expectation and the variance of a random 
variable ^. The diagonal entries 

{WN)^^, 1<1<N, (1.3) 

are i.i.d. random variables, independent from the ofF-diagonal entries, such that 

nWN)^^ = 0, ^{Wn)u = d?, 1<1<N. (1.4) 

We will denote the probability distribution of -^(VFjv)ii by ^i. 

In a similar fashion, in the Hermitian case, we assume that the ofF-diagonal 
entries 

V2dKz{WN)jk,V23m{WN)jk, l<j<k<N, (1.5) 
are i.i.d. centered random variables with probability distribution /i with variance 
and finite fourth moment 7714. The diagonal entries 

{WnU, l<i<N, (1.6) 

are i.i.d. random variables, independent from the off-diagonal entries, with proba- 
bility distribution /ii and finite second moment. 
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While the independence of the matrix entries {WN)ij, ^<i<j<N is crucial 
in our analysis, the requirement that the entries are identically distributed can 
be replaced by certain Lindeberg-Feller type conditions for the fourth moments of 
marginal distributions ([[30]). 

Given a real symmetric (Hermitian) matrix B of order N, we define its empirical 
distribution of the eigenvalues as /is = jfJ2i=i^^iy where Ai < ... < Xn are 
the (ordered) eigenvalues of B. One of the fundamental results of random matrix 
theory is the celebrated Wigner semicircle law (see e.g. [7], [1], [3]). It states that 
almost surely (J-Xn converges weakly to the nonrandom limiting distribution /i^c 
whose density is given by 

d^isc , ^ 1 



(-) = ^v/4^^1[_,.,.,(-). (1-7) 

In other words, for any bounded continuous test function (/3 : M — > K, the linear 
statistic 

1 ^ 1 

-^^(A,) = -Tr(^(XA,)) tr^(^(X^)) 

1=1 

converges to J (p{x)d^sc{dx) almost surely; here and throughout the paper we use 
the notation trjv = "^Tr to denote the normalized trace. 
The Stieltjes transform of the semi-circle law is 



9Az) ■= I = . e C\[-2a,2a]. (1.8) 



It is the solution to 



,2 „2 / 



<j'gi{z)~zg^iz) + l = (1.9) 

that decays to as \z\ — ^ oo. 

In this paper, we are interested in studying the joint distribution of matrix en- 
tries of regular functions of a Wigner random matrix X^. In [28j . Lytova and Pastur 
studied the limit of the one dimensional distribution of y/N {f{XN)ij — E{f{Xiy)ij)) 
in the case of GOE(GUE) ensembles (so the marginal distribution /i of matrix en- 
tries is Gaussian) provided f{x) is a bounded differentiable function with bounded 
derivative. Namely, they prove that 

VN{f{XNh-miXNh))^NiO,^^j^uj\f)), (1.10) 

with /? = 1(2) in the GOE (GUE) case, 

1 Z-^- , 1 



^'(/) := mm = 77 / / (/(^) ~ /(y))'T^ V4a2 - xW^cr^ ' v'dxdy 



l-2aJ-2a ^" 

(1.11) 

where r] is distributed according to the Wigner semicircle law (|1.7p . In the case of 
the off-diagonal entries in the GUE case, the r.h.s. in (I1.10[) should be understood 
as a complex Gaussian distribution with independent identically distributed real 
and imaginary parts, each with the variance ia;^(/). The proof in [2 8) relies on 
the orthogonal (unitary) invariance of the GOE (GUE) ensembles. 

We extend the results of [28] in the following way. We study the joint distribution 
of any finite number of the the matrix entries f{XN)ij The limiting distribution on 
the r.h.s. of ()1.10|) is, in general, no longer Gaussian. Instead, it is the distribution 
of a linear combination of (WN)ij and a Gaussian random variable, independent 
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from {WN)ij (see Theorems 12.31 and 12.71 below). We refer the reader to Remark 
12.41 after Theorem 12.31 for the discussion on when one of the two components in 
the hnear combination vanishes. In particular, the limiting distribution of ^/N x 
{f{X]s[)ij ~ E(/(XAr)y )) is Gaussian if and only if either the marginal distribution 

is Gaussian or /_^2(t — x^dx = 0. 

Our approach requires that / has four continuous derivatives in a neighborhood 
of the support of the Wigner semicircle law, [—2a, 2a]. If the marginal distributions 
H and /ii satisfy a Poincare inequality (j2.15|) then our results hold provided / is 
Lipschitz continuous in a neighborhood of [—2a, 2a]. 

The problem about the fluctuation of the entries of /(Xn) is interesting in its 
own right. However, for us the main motivation to study the problem came from the 
question about the limiting distribution of the outliers of finite rank perturbations 
of standard Wigner matrices (see e.g. [12], [13], ^3\). Let Mn = Xm + Cn, 
where Xn is a random real symmetric (Hermitian) Wigner matrix defined above 
and Cat is a deterministic real symmetric (Hermitian) matrix of finite rank k with 
fixed non-zero eigenvalues Xi(Cn), ■ • ■ , Afc(Cjv) and the corresponding orthonormal 
eigenvectors ui, . . . , Ufc. By the interlacing property, Mn has at most k eigenvalues 
(called outliers) that stay outside the support of the Wigner semicircle law in the 
limit of large N. Capitaine, Donati-Martin, and Feral proved in [13] that the limiting 
distribution of the outliers depends on the localization/delocalization properties of 
the eigenvectors vi, . . . ,Vk. In particular, if the eigenvectors are localized (so only a 
finite number of coordinates are non-zero as iV — >■ oo), then the limiting distribution 
of the outliers is non-universal and depends on the marginal distribution of the 
matrix entries of Wn . 

The results in [13) are proved under the assumption that the marginal distribu- 
tion of the i.i.d. entries of Wn is symmetric and satisfies the Poincare inequality 
(j2.15p . In [33], we have extended the results of [13] to the case of a finite fifth 
moment. Our approach relies on an ideas from [8]. In particular, an important 
step of the proof is the study of the limiting distribution of the eigenvalues of the 
kxk matrix {{vi, RN{z)vj))i<ij<k, where Rn{z) is the resolvent of X^r and (•, •) is 
the standard inner product in C''^ (see Proposition 1 in [33]). Thus, in the localized 
case one is interested in the joint distribution of a finite number of resolvent entries 
of a standard Wigner matrix. 

The rest of this section is devoted to the explanation of the main idea of the 
proof. Complete formulations of the results are given in Section [2j We restrict our 
attention to the real symmetric case since the arguments in the Hermitian case are 
very similar. 

We start by considering the test functions of the form f{x) — which corre- 
sponds to studying the resolvent entries. Define 

Rn{z) := {zIn - Xn)-\ 

the resolvent of a real symmetric Wigner matrix Xm = '^^n for z outside the 
spectrum oi X^. For simplicity, we will consider here a diagonal entry of Rn{z), 
say the (A^, A^)— th entry {Rn{z))nn- The off-diagonal entries can be treated in a 
similar way (see Section [S]) . When it does not lead to ambiguity, we will use the 
shorthand notations Rij for {RN{z))ij and Wij for [WN)ij- Further, assume that z 
is fixed and 3m z ^ 0. By Gramer's rule, the (iV, A^)— th entry of the resolvent can 
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be written as 

Rnn = {z- 4fT^NN - ^b'Gb)-\ (1.12) 



where b is the {N — 1)— dimensional column vector with the coordinates Wn, 1 < 
i < N — 1, 6* its transpose, and G = G{z) is the resolvent of the (iV — 1) x (iV — 1) 
upper- left submatrix X = (^-^Wij)i<ij<N-i of the Wigner matrix Xn, i-e. 

G := {zIn-i-X)-\ 
We note that X can be viewed as a standard {N — 1) x [N — 1) real symmetric 



Wigner matrix since the normahzation by VW instead of \/N — 1 does not make 
any difference in the limit of large N. It is very important in our analysis that 
the random variables Wnn and 5*G5 in (|1.12p are independent. Moreover, in the 
quadratic form 

5*G5= G^jW.iWju (1-13) 

l<jj"<Ar-l 

the vector b is independent from the matrix R. By subtracting and adding E(-^ x 
6*G5) = (T^E(-i-TrG) in the denominator, we rewrite (|1.12l) as 

Rnn ^(z- n^b'Gb) - A^Wnn - i^b'Gb - K{^b'Gb))\ (1.14) 



= [z- a^^TTG) - -^Wnn - (^6*G5 - ^E(TrG)) ) (1.15) 



= [z- a2E(lTri?A.(z)) + 0(1) - -^Wnn - {^b'Gb ~ ^E(TrG)) 



-1 



(1.16) 

where in (|1.16p we used the interlacing property satisfied by the eigenvalues of X^ 
and its submatrix X to write 

E(lTrG) = Ei^TvR^iz)) + O(^) 

since 3m z 7^ is fixed. 

It follows from the semicircle law that 

lim E^Tri?Ar(z) = go-(z), 

A^— foo iV 

where the Stieltjes transform ga{z) of the semicircle law has been defined in (jl.Sp . 
It follows from the calculations in Section [3] (see p.ip in Proposition 13. ip that for 
fixed z 

E^TiR^{z)=gM) + 0{^). 

Therefore, 

Rnn= (z-a^g,{z) + 0{^)-^(wNN + 4^{b*Gb-a^EiTrG)) 

\ -'V -v/iV V V -'V 



(1.17) 

As we have already remarked, the random variables Wnn and (b*Gb — (T-^E(TrG)) 
are independent. The crucial step in the analysis of the fluctuation of the resolvent 
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entries is to prove that 

^ {b'Gb-a^E{TTG)) (1.18) 



converges in distribution as iV — >■ oo to a centralized complex Gaussian random 
variable. We discuss why such a convergence in distribution takes place a few 
paragraphs below but first we note that the Central Limit Theorem for p.lSp 
together with the formula (11.17^ for Rnn and 

z-a g„{z) = — — 

from (|1.9p immediately implies that the normalized resolvent entry ^/N{Rnn — 
ga{z)) converges in distribution as — cx) to the law of a linear combination of 
Wnn and a complex centralized Gaussian random variable. The coefficients of 
the linear combination and the covariance matrix of the complex Gaussian random 
variable are easily computable (see Theorem 12. II in the next section). 

In order to get insight into the limiting distribution of (|1.18l) . it is useful to 
consider first the case of a quadratic form with deterministic coefficients. Let b be 
a random n— dimensional vector with centralized i.i.d. real components with unit 
variance and finite fourth moment, and A„ be an n x n deterministic real symmetric 
matrix such that 

the operator norm ||A„|| < a for all n > 1, (1-19) 

-Tr{Al) 02 as n oo, (1.20) 
n 

1 



n ^ 



y](A„)^j flj as n oo, (1-21) 

=1 

where a > is some constant that does not depend on n. The CLT for x 

(b'^Anb — o'^E(Trv4„)) was first established by Sevast'yanov ([36]) in the case when 
the coordinates of b are i.i.d. standard Gaussian random variables. For subsequent 
developments, we refer the reader to [43], [9], and [5]. In particular, since 

b*Anb ~ a^E{TrA^) = ^ au{bj - ^r") + E ^^^'^J = E ^" 

we can write b*Anb — (T^E(TrA„) as a sum of martingale differences with respect 
to the filtration Tj = cr(6i, . . . , bj), j = 1, . . . , n. It is not difhcult to prove that 
the conditions (|1. 19111. 2T|) imply that the Central Limit Theorem for martingale 
differences (see e.g. [ISj) is applicable and the normalized random variable -^{b* x 
Anb — TrAn) converges in distribution to A^(0, K^al + 02) as n — >■ 00, where K4 is the 
fourth cumulant of the marginal one-dimensional distribution of b. In (|1.18|) . the 
quadratic form is associated with a complex symmetric random matrix G. Thus, 
first of all, one has to study the joint distribution of 

^ (6* «Hc G6 - cr^E(Ti- D\e G)) and ^ (6* 3m G6 - a^E(Ti 3m G)) , 



where 

meG ^ (JReGij, I <i,j < N JmG = {JmGij, I <i,j < N 
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This corresponds to the choice of 

An = xmeG + y3mG, N (1.22) 

where x and y are arbitrary real numbers. The second difference is that An in ()1.22|) 
is random. However, the CLT still holds if ||A„|| is bounded from above by a non- 
random constant a with probability 1 and iTr(A^) and ;^ X]r=i{^")?i converge 
in probability to non-random limits (see e.g. ^ and the Appendix by J.Baik 
and J.Silverstein in [H]). We note that ||G|| < | . The desired convergence 
in probability will easily follow from the self-averaging property of the resolvent 
entries established in Section [3] (see Proposition 13. 1|) . The generalization of the 
CLT to random An is not unexpected since the distribution of -^{b*'Anb — TrA„) 

conditioned on the matrix entries of An is approximately -/V(0, ^"^j^(A„)fj -|- 
iTr(A^)) for large n, and the expression in the variance converges to a non-random 
limit as n 00. 

As a result, one obtains that the term {UGh — (T^E(TrG')) in (|1.17p converges 
in distribution as 00 to a complex Gaussian random variable with zero math- 
ematical expectation and the covariance matrix given by the r.h.s. of (|2.7ll2.9p with 

w — z and Lp^^[z,w), Lp {z,'w)^ and ipj^ {z,w) defined in (j2.2ll2.5|) in Section[5] 

Since it is independent from Wnn^ the limiting distribution of \/N(Rnj\i — ga{z)) 
is given go-(z)^ multiplied by the convolution of the marginal diagonal distribution 
fj,i and the complex Gaussian. 

To study the joint distribution of several resolvent entries Ri^j^{zi), 1 < I < k, 
one follows a similar route. The main new ingredient is a multi-dimensional GLT 
for random bilinear (sesquilinear in the complex case) forms (see Theorem 6.4. in 
[8] and Theorem 7.3 in for the convenience of the reader we reproduce the last 
one as Theorem 17. II in Section [T]). Thus, one is able to prove the result of Theorem 
12.31 in Section [5] for the test functions of the form 
k ^ 

/(x)= Vcfc , CfeGR, ZkeC, Jmzk^O, 1 < z < fc. (1.23) 

^ Zk-x 

The finite fourth moment condition on the off-diagonal entries of Wn and the 
finite second moment condition for the diagonal entries of Wn imply that HA^Arll — >■ 
2a a.s. a.s N —i' 00 ( 6 ). Therefore, the limiting fluctuations of normalized matrix 
entries of /{Xn) do not change if we alter / outside [—2a — S,2a + 6], where <5 > 
is an arbitrary fixed positive number. In particular, one can replace / by fh where 
h e C°°(M) is a function with compact support such that h — 1 on [—AI,M] for 
sufficiently large M. 

In order to extend the result of Theorem 12.31 from (|1.23p to more general test 
functions /, one approximates / by test functions /,„ of the form (jl.23p so that 
11/ ~ /mil — >■ as m — >■ 00 in an appropriate norm and NY{f{Xpf)ij — fm{XN)ij) 
goes to zero uniformly in N when m — )> cxd, where V denotes the variance. 

This program is the easiest to implement when the marginal distributions of the 
entries of Wn satisfy the Poincare inequality. Indeed, for a Lipschitz test function 
/, a matrix entry f{X)ij is a Lipschitz function of the matrix entries of X (see e.g. 
[14j). Therefore, as a direct consequence of the Poincare inequality for the marginal 
distributions of Wn one gets the bound 

nfi^Nh) < (1.24) 
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where |/|£ sup^^y 




is the Lipschitz constant. We note that N x 



W{f{XN)ij) goes to zero if |/|£ goes to zero. Approximating a Lipschitz continuous 
compactly supported test function / by functions fmh, where fm, m > 1, are of 
the form (|1.23p one finishes the proof. 

If the marginal distributions do not satisfy the Poincare inequality, one needs to 
impose some additional smoothness condition on / to obtain an analogue of (jl.24p . 
We refer the reader to the bound (|2.22p in Proposition 12 . II in Section [2] The proof 
of (|2.22p consists of two steps. First, one estimates the variance of the resolvent 
entries and proves 



where Pg is some polynomial of degree 6 with fixed positive coefficients (see (13. 3p 
in Proposition 13.11 in Section (H]). The proof of this bound is a bit long but quite 
standard and relies on the resolvent technique (see e.g. [3], [H], [35], [37], and [37]). 
In particular, many computations are similar to those used in the derivation of the 
master loop equation in the proof of the Wigner semicircle law by the resolvent 
method. 

To extend the last estimate to more general test functions we use the Helffer- 
Sjostrand functional calculus discussed in Section S] (see e.g. [H], [T^, [3^, or the 
proof of Lemma 5.5.5 in [T]). One then requires that / have four derivatives to 
compensate for the | Umz|~^ factor in the upper bound on Vi?^ (z). 

The rest of the paper is organized as follows. Section [2] is devoted to formulation 
of our main results. Perhaps, during the first glance at the paper the reader could 
just look at Theorem [23] in the real symmetric case (the analogue in the Hermitian 
case is Theorem 12. 7|) and omit the rest of the section. Theorems 12.11 and 12.51 
deal with the resolvent entries and are important building blocks in the proofs of 
Theorems 113] and I2I7J Theorems |12] HI] and |22] prove slightly stronger results 
since as they assume that the marginal distributions satisfy the Poincare inequality. 

The actual proof starts in Section |3] which is devoted to estimates on the math- 
ematical expectation and the variance of the resolvent entries. The main results 
of this section are collected in Proposition 13.11 During the first reading of the pa- 
per, the reader might wish to skip long but rather straightforward computations 
in Section [3] and jump to the next section once the statement of Proposition 13 . 1 1 is 
absorbed. 

In Section |4] we extend our estimates to the matrix entries f{XN)ij for suffi- 
ciently nice test functions / by applying the Helffer-Sjostrand functional calculus. 
The main result of Section 2] is the proof of Proposition 12.11 

Section [5] is devoted to studying the fluctuation of the resolvent entries and 
contains the proofs of Theorems 12. 1[ 12. 2[ 12.51 and 12.61 The proofs of Theorems 12.11 
and 12.61 follow the route explained above when we discussed the fluctuation of a 
diagonal resolvent entry. The proofs of Theorems 12.21 and 12.61 could be omitted at 
the first reading as they prove the functional convergence in a special case when 
the marginal distributions satisfy the Poincare inequality. 

Theorems 12.31 and 12.41 are proved in Section |6] The proofs of the corresponding 
results in the Hermitian case, namely Theorems 12.71 and 12.81 is very similar and 
mostly left to the reader. 

In the Appendix, we discuss various tools used throughout the paper. 
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We will denote throughout the paper by consti,Consti, various positive con- 
stants that may change from line to line. Occasionally, we will drop the dependence 
on TV in the notations for the matrix entries. 

We would like to thank M. Shcherbina for useful discussions and S. O'Rourke 
for careful reading of the manuscript and useful remarks. 



First, we consider the resolvent entries. In Theorems 12.11 and 12.51 formulated 
below, we study the limiting joint distribution of a finite number of resolvent matrix 
entries of a real symmetric (Hcrmitian) Wigner matrix Xjs;. We recall that the 
resolvent Rn{z) of X^^ is defined as 



for z outside the spectrum of Xj^. We will be interested in the limiting joint distri- 
bution of a finite number of the resolvent entries. Since the entries of Xj^ are i.i.d. 
random variables up from the diagonal, we can study, without a loss of generality, 
the joint distribution of the resolvent entries in an to x m upper-left corner of the 
matrix provided that m is an arbitrary fixed positive integer. 

Let us denote by R^"^\z) the m x m upper- left corner of the matrix Rn{z). 
In a similar fashion, we denote by W'^"'\ X^™), the m X m upper-left corner of 
matrices Wn and Xn, respectively. If the reader is put off by some cumbersome 
formulas/notations in this section, he/she can always assume to = 1 and deal with 
just one diagonal resolvent entry. The case to > 1 does not require any significant 
new ideas. 

Consider a matrix- valued random field 



Clearly, Tjv(z) is a random function on C \ [— 2cr, 2(t] with values in the space of 
complex symmetric to x to matrices {Tn{x) is real symmetric for real x). Let us 
define 



2. Formulation of Main Results 



Rn{z) := {zIn-Xn)-\ 




(2.1) 



1 1 1 



■\/ 4(T^ — x^dx 



if W z, 
if w = z. 



(f{z,w) : 



z ~ X w ~ X 27r(T^ 




(2.2) 
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for z,w e C \ [—2a, 2a]. Clearly, ip{z,w) — E , where 77 is distributed 



according to the Wigner semicircle law (|1.7p . We also define 

Lp++(z,w) := [ Die 9\z y/ia"^ - x^dx (2.3) 

J_2a z — X w — X zvrcr^ 



= ^ ('^(^: W) + ip{z, W) + (^(2;, W) + lf{z, W)) , 

3m 3m ^V^tr^T^da; (2.4) 

_2cr z — X w ~ X 27rcr^ 

= (</'(^j w) + fiz, w) - (p{z, w) ~ <^(z, w)) , 

/2a 2^ 2^ 
9le 3m — :T\/4a^ - x^dx (2.5) 
_2ct z — x w — X 2'Ka'^ 

= (95(2:, ■*«) - ¥'(5,«)) - '^{z,w)) , 

Theorem 2.1. Let = "^^'^^ random real symmetric Wigner matrix hl.lV 
QT^p. T/ie random field Tn{z) in h2.1\) converges in finite- dimensional distributions 
to a random field 

T(z)^.92(z)(Ty(™)+y(z)), (2.6) 

where l^^™) is the mx m upper-left corner submatrix of a Wigner matrix Wn, and 
Y{z) = {Yij{z)) ,Yij{z) = Yji{z), 1 < i,j < m, is a Gaussian random field such 
that 

CovCO\eYjj{z),meYjj{w)) = K4{p)mc ga{z)me g^{w) 2a'^ip++{z,w), (2.7) 

Cov {3m Yjj{z), 3m Yjj{w)) — n^^fj,) 3m ga{z) 3m ga{w) +2a'^ if {z,w), (2.8) 

Cov{yieYjj{z),3mYjj{w)) = K4{n)meg^{z) 3m g„{w) + 2a'^ip+^{z,w), (2.9) 
CovimeY,,{z),meY,,{w))^a^ip++{z,w), i^j, (2.10) 
Cov{3mY,j{z),3mY,j{wj) ^a'^Lp^^{z,w), i^j, (2.11) 
<Cov{^zY^j(z),3mY,j(w)) = a^Lp+^(z,w), i ^ j, (2.12) 

where the fourth cumulant K4(/i) := / x'^fj,{dx) — 3(/ x'^ fi{dx))'^ ~ 7714 — 3cr^. 

In addition, for any finite r > 1, the entries Yi^j^ (zi), 1 < ii < ji < m, I < I < r, 
are independent provided (iinjii) 7^ {ii2iji2) f'^''" ^ 1^ h ^ h l£ r. 

Remark 2.1. If z — x ^ M. \ [—2a, 2a], then Yij{x), 1 < i < j < rn, are 

independent centered real Gaussian random variables with the variance given by 

Y{Yu{x)) = n^{ii)gl{x) - 2a^g',{x), I < i < m, (2.13) 
Y{Y,,{x)) = -a^g',{x), l<i<j<m. (2.14) 

Let I? C C \ [—2a, 2a] be a compact set. The distribution of Tjv(z),z G T>, 
defines a probability measure V^f on C(2?, C™^™"''^-'/^). One can prove functional 
convergence in distribution for the random field Tjf(z), z £ T>, provided /i and 
fii satisfy some additional conditions on the decay of their tail distributions. For 
simplicity, we will consider the case when and /xi satisfy the Poincare inequality. 
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We recall that a probability measure P on M*^ satisfies the Poincare inequality 
with constant u > if, for all continuously differentiable functions / : M*^ — > C, 

Vp(/) = Ep - Ep(/(x))|2) < iEp[|V/(.T)|2] (2.15) 

Note that the Poincare inequality tensorizes and the probability measures satisfying 
the Poincare inequality have sub-exponential tails ([H], [T]) . By a standard scaling 
argument, we note that if the marginal distributions fi and ni of the matrix entries 
of Wn satisfy the Poincare inequality with constant v > then the marginal 
distributions of the matrix entries of = '^^n satisfy the Poincare inequality 
with constant Nv. 

Theorem 2.2. Let = "^^a' a random real symmetric Wigner matrix 
and the marginal distributions jj, and fii satisfy the Poincare inequality 
S2.15]) . The probability measure Vj^ on (7(2?, C"^"*"*"^-*^^) given by the random field 
Tn{z) in 12.1]) weakly converges to the distribution of the random field T (z) defined 
in Theorem ] 2. 1\ 

Next, we extend the results of Theorem 12. ll to the matrix entries of /(X^r) for 
regular /. We say that a function / : M — )• M belongs to C" (/) , if / and its first n 
derivatives are continuous on /. We will use the notation C"(M) for the space of n 
times continuously differentiable functions on R with compact support. We define 
the norm on C"{I) for compact / C M as 

||/||c"(/) := max(|dV/da:'=(a:)|, 0<k<n, xel). (2.16) 
We also define for / e C"(M) 

/CO 
\d''f/dx''(x)\dx, 0<k<n), (2.17) 
-oo 

and 

||/||„4,+ := max (^J^{\x\ + l)\^{x)\dx, < ? < . (2.18) 

Clearly, the right hand sides of (|2.17p and (|2.18p could be infinite. 

We start with Proposition 12.11 that holds both in the real symmetric and Her- 
mitian cases. 

Proposition 2.1. Let — '^^n be a random real symmetric (Hermitian) 
Wigner matrix defined in hLl\[l.Ji^ i^ jj. 51177^) ). 

(i) Let L he some positive number and f S Cj(R) with supp{f) C [— 
Then there exists a constant Const{L, /j,, ni) that depends on and /ii, such 

that 

\E{f{XNh)- r f{x)-^^A<T^-xHx\ < Const(L,/.,/ii)^^^^itMl), i < ^ < 



(ii) Let f e C^(M), then 



-2a 
8/ 



(2.19) 



\&{f{XN)^^)- / f{x)^^^a-^~xUx\ (2.20) 



-2cr 



< Const(/^, ^i)Mli|L±^ l<i<N. 
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YifiXNh) < Consti^^,^l,y-^, 1 < i,j < N. (2.22) 



(Hi) Let f £ C^(M), then there exists a constant Const{^, fix) such that 

|E(/(Xjv),fc)| < Const{^l,^l,)Mhl^ l<j<k<N. (2.21) 

(iv) Let f £ C"*(R), then there exists a constant Const{p, such that 

1^1 
N 

(v) Lf fj, has finite fifth moment, /ii has finite third moment, and f £ C^'^{R), 
then one can improve 112. 2 namely 

|E(/(Xw),fc)| < Const{p,^ll)^^, 1<J <k<N. (2.23) 

Remark 2.2. In f33', we extend the results of Propositions \2.1\ to the case of 
{u'^^y f {X j\i)v'^^^) , where w^^^ and ii^^^ are arbitrary nonrandom vectors from 

In order to formulate our next theorem, we need to introduce several notations. 
Recall that we defined w^(/) in (II. lip as 

^'(/) = V(/(77)) = i r r {fix) - f{y)f-^JAa^ - xV4ct2 - V^dxdy, 



^ J -2a J -2a " ' " 

where rj is distributed according to the Wigner semicircle law (jl.7p . In addition, 
we define 

a{f) E (f{r^)'A = - T xf{x)-^^Aa^-xUx (2.24) 

and 

/?(/) (^/(,7)^^^^) - ^ f_J{x){x^ - (2.25) 

We recall that a function / : M ^ M is called Lipschitz continuous on an interval 
/ C M if there exists a constant L such that 

\f{x)- f{y)\<L\x-yl for all ye/. (2.26) 

We define 

|/k« = sup^^M^, (2.27) 

x^y \x-y\ 

and 

Ifl .„T, \f{x)- fiy)\ .r,<,oN 

\f\c,5= sup j j . (2.28) 

x^y., x,ye[-2<y-&,2a+S] F — y| 

Finally, let us introduce a C°°(M) function h{x) with compact support such that 

h{x) = 1 for X £ [-2(7 -5,2a + 5], 5>Q. (2.29) 

Theorem 2.3. Let Xn = --j=Wm he a random real symmetric Wigner matrix \L1V 
I J./^[ ). Let / : R — > M he four times continuously differentiahle on [—2cr — (5, 2(7 + 8\ 
for some 5 > Q. Then the following holds, 
(i) For i = j, 

VN{f{XN)^^-^{{fh){XN)^.i)) (2.30) 
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converges in distribution to the sum of two independent random variables 
^^^Wii and N(0,vf{f)), where h is an arbitrary C^(R) function satisfying h2.29]) . 
and 

vlif) 2 ^^Xf) - a'if) + ■ (2.31) 

/// is seven times continuously difjerentiable on [—2(7— S, 2a +S], then the statement 
still holds if one replaces M {{fh){XN)ii) in h2. 30]) by 

f.2cr -, 

/Wt; — :tv/4(t2 -a;2da;. (2.32) 
2a 27ra2 

(ii) For i ^ j, 

VN(f{XNh-E{{fh){XNh)) (2.33) 
converges in distribution to the sum of two independent random variables 
^W,j and N{0,d^{f)), with 

d\f):=Lo\f)-a'[f). (2.34) 

/// is six times continuously differentiable on [—2a — d, 2a + d], then one can replace 
E{{fh){XNhj) m ^EM by 0. 

(Hi) For any finite m, the normalized matrix entries 



'N ifiX^h - mfh){XNh)) , 1 < z < J < m, (2.35) 
are independent in the limit N ~> oo. 

We follow Theorem 12.31 with several remarks. 



Remark 2.3. If f e C^(R), and \\f\\4,i < oo, where \\f\\4,i is defined in |g.J7| ), 
then it follows from Proposition \2.1\ that the centralizing constants in \2.3U\) and 
V2.33\) can be taken to be E{f{XN))^j. 

Remark 2.4. It follows from the definition of the fourth cumulant that "^^4-* > —1, 

2 2 

with the equality taking place only when fi is Bernoulli. Since 1, ^, and ^ ~'f , 
are the first three orthonormal polynomials associated with the semicircle measure 
Jj. 7| ), it follows from the Bessel inequality that for the diagonal entries the variance 
v\{f) of the Gaussian component in \2. 31\) is zero if and only if f(x) is either a 
linear function of x or f{x) a quadratic polynomial and fi Bernoulli. Similarly, for 
the off-diagonal entries one has that the variance d^{f) of the Gaussian component 
in \2. 34\ ) is zero if and only if f(x) is linear. 

The statement of the Theorem \2.3\ also implies that the limiting distribution of 
the normalized {ij)th entry of f(XN) is Gaussian if and only if either fi for i ^ j 
(correspondingly fj,i for i = j) is Gaussian or ¥,(7] f (r])) = 0. The same holds in the 
Hermitian case. 

Remark 2.5. Utilizing Proposition 1 in [37] . one can extend the result of Theorem 
\2.3\ to the test functions satisfying 

11/11' = /(I + \k\f'\fm^dk < 00, 5 > 3, />) := J e-^^^^f{x)dx, (2.36) 

and, more generally, to the functions that coincide with the functions in 112. 36]) on 
[-2a -d,2a + S] for some S>0 (^). 
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Remark 2.6. There has been a significant body of work on the Central Limit 
Theorem for TrfiXn) = E^i fi^N)u- We refer the reader to [2], [5], [4], [23, 
and [37 . and the references therein. In particular, in [37] . the CLT is proved 
assuming that the fourth moment of the marginal distribution is finite and \\f\\s < oo 
for s > 3/2. 

If /i and 111 satisfy the Poincare inequality, one can prove convergence in distri- 
bution for the matrix entries of Lipschitz continuous test functions. 

Theorem 2.4. Let Xn — '^^n be a random real symmetric Wigner matrix 



and the marginal distributions fi and /ii satisfy the Poincare inequality 
f2.15\) . Then the following holds. 

(i) If f : 'M. ^ R is Lipschitz continuous on [—2(7 — (5, 2(t + 5] that satisfies a 
sub- exponential growth condition 

< aexp(6|x|) for all x £ R (2.37) 

for some positive constants a and b, then the results of Theorem \2.3\ hold with the 




centralizing constants K(f{XN)ij)) in S2. 30\) and S2.33\) . 
Moreover, 

P{\f{XN)^J~E{f{XNh)\ > t) (2.38) 

/ \/vNt \ 
< 2K exp I - \+{2K + o(l)) exp 

where \f\c,S is defined in i2.28\) . 

A' = -^2Mog(l-2-i4-*), (2.39) 

i>0 

and V is the constant in the Poincare inequality 12.15\) . 

(a) If f e C^(R) (correspondingly, f e C^(R), / £ C^"iR)) and f satisfies the 
subexponential growth condition {2.31^ , then the estimate \2.19\) (correspondingly, 
ifOTl) . IMM ) from Proposition [O holds. 

(Hi) If f is a Lipschitz continuous function on M, then 

P{\f{XNh--m{XNh)\ > t) (2.40) 



< 2ifexp 



'vNt 



where \f\cM is defined in 1^2.27^ . 

Remark 2.7. If f is a Lipschitz continuous function on [—2a — 6,2a + S] that does 
not satisfy the subexponential growth condition ^2.37\ l then the results of Theorem 
'YMstill hold with the centralizing constants E {{fh){XN)ij))) in 112.30]) and 112.33]) . 

In the second part of this section, we formulate the analogous results in the the 
Hermitian case. 

As in the real symmetric case, consider a matrix- valued random field 



- N 



(z) - Vn (i?(")(z) - g,(z)/„) , zeC\ [-2(7, 2a] 



Clearly, Tn{z) is a random function on C \ [—2a, 2a] with values in the space 
of complex m x m matrices. Tat (a;) is Hermitian for real x and, more generally, 

TAr(z) = TAr(z)*. 
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Theorem 2.5. Let Xn — -^Wn be a Hermitian Wigner matrix il.5}\1.6]) . The 

random field Tjv(z) converges in finite- dimensional distributions to a random field 

T{z)^gl{z){W^"^^+Y{z)), (2.41) 

where is the mx m upper-left; corner submatrix of a Wigner matrix Wm, and 

Y{z) = {Yij{z)) , 1 < j, J < TO, is a Gaussian random field such that 



Cov{mtY,,iz[ 


),9{tY,,{w)) = 


Kiip.) d\cga{z) yieg^iw) - 






w), 


(2.42) 


Cov (am Yjj {z),3mYjj{w)) ^ 


- Ki{^) 3m ga{z) 'Imgai'w) 


+ aV 


-{^ 




(2.43) 


CovC^RcYjjiz 


),3mYjj{w)) = 


K^ifj.) Ulega{z) Jm ga{w) 


+ aV- 




w), 


(2.44) 


Cov{meY,j{z] 


\,mcY,,{w)) = 


^<^'^{'P++{z,w) + ip—{z, 


w)), i 






(2.45) 


Cov{3mY^j{z 


),3mKy(u))) = 


^<^'^{f++{z,w) + if—{z, 


w)), i 






(2.46) 


Cov{meY,j{z] 


\,3mY^j{w)) = 


icr*(^+_(z,-U7) - Lp+^{w 


,z)), i 


+ 3- 




(2.47) 



where the fourth cumulant K4{fi) :— E,\{Wn)i2\'*' ~ '2E,\{Wn)i2\'^ — — 2a^. 

In addition, for any finite r > 1, the entries Y^j^ (zi), \ < ii < ji < m, 1 < I < r, 
are independent provided (iinjii) 7^ {ih^ih) /"^ 1 < 7^ Z2 < 



Remark 2.8. If z = x e R\[-2(7, 2cr], thenYa{x), I < I < m, D\cYij{x), JmY.jix), 
^ l£ i < j l£ "Tti are independent centered real Gaussian random variables with the 
variance given by 

YiYuix)) = Ki{fi)glix) - a^g'Jx), 1 < / < m, (2.48) 

YimcY^.ix)) = -^a^g'A^), ¥{3mY,,ix)) = -^a^g'J^), 1 < i < j < m, (2.49) 
Cov{D\tY^j{x),3mY,j{x)) ==0, 1 < i < j < m. (2.50) 

Let V, as before, be a compact subset of C \ [—2a-, 2a-]. The distribution of 
TN{z),ze V, defines a probability measure Vj~^ on C(X', C"(™+i)/2). 

Theorem 2.6. Let Xn — '^^^n be a random Hermitian Wigner matrix lil.5\ - 
\L6\) and the marginal distributions fj, and /xi satisfy the Poincare inequality i2.15\) . 
The probability measure Vj^f on (7(2?, C™*-™"*"^^/^) given by the random field Tn{z) 
weakly converges to the distribution of the random field T(z) defined in Theorem 
\M 

Next theorem extends the results of Theorems 12 .31 and 12 .41 to the Hermitian case. 
We recall that uj^{f),a(f), (i{f), and d'^{f) have been defined in (fLTTll . ([221, 
(12251), and (12:3411 . 

Theorem 2.7. Let Xj^i — -^Wjy be a random Hermitian Wigner matrix U.51 
\1.6]} . Let / : R — > M be four times continuously differentiable on [—2(7 — (5, 2(7 + S\ 
for some 5 > Q. Then the following holds, 
(i) For i = j, 

VNif{XN)u-Ei{fh){XNh)) (2.51) 
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converges in distribution to the sum of two independent random variables 
^^^Wii and N{0,V2{f)), where h is an arbitrary C^(R) function satisfying h2.29\) . 

4i^f):=^\f)-a\f) + ^p\f), (2.52) 

and 

h4p) - E|(W^jv)i2|^ - 2E|(M^Ar)i2p ^ nn - 2a\ 

If f is seven times continuously differ entiable on [—2a — S,2a + d], then one 
replace E{{fh){XNh) in KM) by ^MB- 
(a) For i 7^ J, 



ifiX^h - E {ifh){XMh)) (2.53) 
converges in distribution to the sum of two independent random variables 
^^^Wij and complex Gaussian N(0,d'^{f)) with i.i.d real and imaginary parts 
iV(0, ^c?^(/)). /// is six times continuously differentiable on [—2a — S,2a + S], then 
one can replace E ((//i)(Xjv)y ) in h2. 5S\) by 0. 

(Hi) For any finite m, the normalized matrix entries 



'N ifiXNh - nifh)iXNh)) , l<i<j<m, (2.54) 
are independent in the limit — > oo. 

Theorem 2.8. Let Xjv — '^Jn^^ random Hermitian Wigner matrix il.5l 

\1.6]) and the marginal distributions fi and /ii satisfy the Poincare inequality \2.15]) . 
Then the following holds 

(i) V f is a Lipschitz continuous function on [—2a — 6,2a + 6] that satisfies 
the subexponential growth condition {2.31^ , then the results of {2.1^ hold with the 
centralizing constants K {{fh){XN)ij) in i2.51\) and i2.53\) . 

Moreover, 

F{\f{XNh,-m{XNh)\ > t) (2.55) 

where K is as in i2.39\) . v is the constant in the Poincare inequality h2.15\) . and 
\f\c.S defined m [tM) . 

(ii) If f € C^(K) (correspondingly, f € C^{R), f e Ci°(M)) and f satisfies the 
subexponential growth condition {2.37) , then the estimate i2.19\) (correspondingly, 
Ji2Jl\) . K2^) } from Proposition\2^ holds. 

(Hi) If the marginal distributions /i and /ii satisfy the Poincare inequality h2. 1 5\) 
and f is a Lipschitz continuous function on M, then 

¥{\f{XN)^,-n.f{XNh)\>t) (2.56) 



where \f\cM is defined in \2.2T^ . 

Most of the proofs will be given in the real symmetric case. The proofs in the 
Hermitian are essentially the same. 
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Remark 2.9. Theorems \2.1l \2.S\ \2.5\ \2. 7| can be extended to the case when the 
matrix entries (W'Af)y, 1 < « < j < o,re independent but not identically dis- 
tributed [30] . In the real symmetric case, one requires that the off-diagonal entries 
satisfy 

E(Wjv),j = 0, ^{{WN)^3) = l<i<j<N, sup E((W^w),y)4 = m4 < oo, 

mi{i) :— Ym\ — E|(WAr)ij|^ exists for I < i < m, 

Ljv(e) — > 0, as N ^ oo, Ve > 0, where 

^^(^) = 1^ E ^{\iWNhW{\{WNh\>eVN)), 

l<i<j<N 

Li,N{() —^0, as N ^ oo, Ve > 0, 1 < i < m, where 
and the diagonal entries satisfy 

nWN)^^ = 0, SUp^,NE\{WN)^^\^ < OO, 

In{^) — > 0, as N oo, Ve > 0, where 

'^W = ^ E ^{\iWNhWmNh\>eVN)). 

l<i<N 

In the Hermitian case, one requires that, in addition, U\e{WN)ij is independent 
from 3xn{WN)ij, ^ < i < j < N, and 

2 

Y{mc(WNh) = Y{JmiWNh) = y , l<t<j<N. 



3. Mathematical Expectation and Variance of Resolvent Entries 

In this section, we estimate matliematical expectation and variance of resolvent 
entries Rij (z) :— {Rpf{z))ij. Wittiout loss of generality, we can restrict our attention 
to the real symmetric case. The proofs in the Hermitian case are very similar. 
Usually, we will assume in our calculations that /ii = fi. The proofs in the case 
fii fj, are very similar. From time to time, we will point out the (small) changes 
in the proofs one needs to make if fii ^ ii. 

Proposition 3.1. Let = '^^n be a random, real symmetric (Hermitian) 
Wigner matrix defined in m.l^l.4\ ) f il.5}\1.6\) ) and Riq{z) = (zIn — Xj^)^^ where 
z € C. We will denote by Pi{x), I > 1, a polynomial of degree I with fixed positive 
coefficients. 
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Then 

EtmRN = ERn{z) = ga{z) + O (j^^^^^^ > (3-1) 

uniformly on bounded subsets o/ C \ M, 
jitn 

V%(z) = O I •^'^'■^^"^^^ — ^ ) , l<i,j <N, uniformly onC\R. (3.3) 



Rij{z) = O ( ^^^^•'""^/^ — ^ ) , 1 < i 7^ j < iV, uniformly on C \ R, (3.2) 



/n addition, if fi has finite fifth moment and fii has finite third moment, then 

ER,j{z) = O (^^-^^^^^ — , l<i^j<N, uniformly onC\R. (3.4) 

Remark 3.1. We refer the reader to [18| (see e.g. Theorem 2.1 there) for the 
optimal (up to logiV factors) estimates on the resolvent entries with the correct 
|3Tnz|~^ behavior. The authors in |18j require that the marginal distributions are 
subexponential. For related results, we refer the reader to the survey [T7|. 

Proof. We use the resolvent identity (|7.2I) to write 

zER,j{z) = d,j +"^E(X^kRkj)- (3.5) 

k 

Applying the decoupling formula (|7.1I) to the term K{XikRkj(z)) in (|3.5p . we obtain 
the equation 

2 

zEi?y(z) = 6,jWE{R^jiz)tTNRNiz)) +—E{{RNiz)%)^— (E(i?,,(z)i?y- (z)))+rA,, 

(3.6) 

where rjv contains the third cumulant term corresponding to p = 2 in (|7.ip for 
k ^ i, and the error terms due to the truncation of the decoupling formula (|7.ip at 
p = 2 for A: 7^ i and at p = for k — i. We rewrite (|3.6p as 

zER,j{z) = S^J + a^ER^j{z)g„{z) + a^ER^j{z) (EtrAri?Ar(z) - g„{z)) (3.7) 
+ a^CoviR,,{z),trNRNiz)) + ^E {{Rn{z)%j) ~ ^ (E(i?,,(z)i?,,(z))) + r^. 



(3.8) 



We claim the following estimates uniformly on C \ . 
Lemma 3.1. 



Co«(i?,,(z), tr^R^{z)) = O (^M^A^ , (3.9) 
^E((i?.(.)%.)=o(^^^^), (3.10) 
..^of^^MlV (3.11) 



V ^ 

Proof. The bound p. 101) immediately follows from (|7.9p . To obtain p.9p . we again 
use (|7.9I) and the estimate 

V(tr^i^^.(z)) = O ( 3^^' ) , (3.12) 
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from Proposition 2 of [37] . It follows from the proof that the bound is valid provided 
the second moments of the diagonal entries are uniformly bounded and the fourth 
moments of the off-diagonal entries are also uniformly bounded ( |38| ) . 

Now, we turn our attention to (|3.1ip . First, we note that the term x 
(E(i?ii(z)i?ij(z))) in p.Sp is O ^ | 2m\[^N ^ which immediately follows from (|7.9p . 

The third cumulant term gives 



-[4E( J2 R^lR^kRkk) + 2E( RuRkkRk,) (3.13) 



2!iV3/2 ' 



2EiY,iR^k)^RJk)]■ 



Since 



k:k^i 

1 1 



E l^^fcl' < ll^ll' < Jj^^ and \R,,\iz) < (3.14) 

we conclude that the third cumulant term contributes O (^ivTlmTp) to rjv in p.Sp . 
In a similar way, the error term that appears due to the truncation of the decoupling 
formula (|7.ip at p = 2 is O (ivyaS^) • In.dced, it can be written as a sum of 0{N) 
terms, where each term is bounded by O {-^ \ 3m zj"**) . Lemma |3. II is proven. □ 
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Proof of / TO) 

Now, we turn our attention to p.ip . For 

gN{z) := EtTNR = Ei?ii 
one can write the Master Equation as 

zgN{z) - l+a^g%{z)+a^Cov{Rn{z),tTNRN{z)) + ^E{{RN{zf}ii)+rN, (3.15) 

by applying (|7.1I) to E,{XikRki{z)) and summing over 1 < fc < iV. As before, 
contains the third cumulant term corresponding to p = 2 in (|7.ip for k i, and the 
error due to the truncation of the decoupling formula (|7.ip at p = 2 for k i and 
at p = for /c = i. 

By p.9|) and (|3.10p . we bound the third and the fourth terms on the r.h.s. of 
(I3T5)) by O (^ ^3(|Jmzr^) ^ Q ^ P2{\0mz\-^) -^ ^ respectively Thus, we obtain 

zg^iz) = 1 + a'g%iz) + O (^M^plH^ , (3.I6) 

uniformly in z G C \ M. 

We now show that the bound (|3.16p implies (j3.ip uniformly in z satisfying 

\z\ < T, and Jmz^O, (3.17) 

where T is an arbitrary large fixed positive number. Our proof follows closely argu- 
ments from [12] Proposition 4.2, ^Tj Section 3.4, and ^20j Lemma 5.5, Proposition 
5.6, and Theorem 5.7. Define 

Qn ^{z:\z\<T + l, I 3mz\ > LN-^'"^}, (3.18) 
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where L > will be chosen to be sufficiently large. Then for z G Qn, the error 
term N\ymz\i = 0(iV-i/^). Therefore 

zgNiz) - a^gUz) = 1 + OiN-'^'), (3.19) 

so 

\gNiz)\ > 5i > on Qn. (3.20) 

Let 

sn{z) = (j'^gNiz) H 

gN[z) 

Then it follows from and (IX^ 

SNiz) z = a'g^iz) + ^ " - = ^(^^3^^^ ^^'^^^ 



on Qat. Since 



cr^57v(2) H ^ a^ga{sNiz)) + 



9n{z) ga{sN{z)) 
we conclude that 

gM{z) = g„{sN{z)), (3.22) 

first for | 3m z \ > v^cr, and then for all z € Qat by the principle of uniqueness of 
analytic continuation. 

Choosing L in p.lSp sufficiently large, we have that 

|3msAr(z)| > ^|3mz| 
on Qa,. Since |%^| < p^^, we conclude that and ((X^ imply on 

If I 3m z| < then ^vfoi;^ ^ ^nd 



3mz| V^l^^nzl' 

Therefore, the estimate p.ip is proven. 
Proo/ o/ iTO) 



Now, we prove p.2p . It follows from p.7ll3.8p and Lemma [3. II that 

^E((P.(.))^). + 0(^ 



zEi?i2 = a2ff^(z)Ei?i2 + '^n{RN{z)f )i2 + O f ^"^'^""/' ) . (3.23) 



Therefore, 



(z - a2.g^(z)) Ei?i2 = O ('^^i^^^^') . (3.24) 



It follows from p.l6p that 

gN{z){z - a^g^iz)) = l + ^^2:^0^^ . (3.25) 
It follows from (1X24)) and (1X25)) that 

('l + 0(^ifl^1V«..-0('M2E£n, (3.26) 



N / / V iV 
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Consider On = {z '■\ 3m z\ > LN ^/*}, where the constant L is chosen sufRciently 
large so that the O (^El!A2l^ — )^ term on the l.h.s. of (|3.26p is at most 1/2 in 
absolute value. Since |f/Ar(z)| < | , we obtain 



,3mz| V N J \ N 
for z G Oat. On the other hand, if | 3m z\ < LN^^^^, then ymz\^ ^ L~^, and 

|Ei?i2| < ttA-t = o( ,J . (3.28) 



Combining p.27p and p.28p . we obtain 
Prooj of iTOI) 

Now we proceed to prove the variance bound p.3p . We apply the resolvent 
identity (|7.2p to write 

zE{R,^{z)R,j{z)) ^ER,j{z)S^j + Y,^iX^kRkJR^J{z)) . (3.29) 



Applying the decoupling formula (|7.ip to the term E {XikRkj{z)Rij{z)) in p.29p . 
we obtain 

zE {R,j{z)R,j{z)) = ER,j{z)S,j + a^E {R,j{z)tTNRN{z)R^jiz)) (3.30) 
+ ^E{{Rn{z)%,R.,{z)) + ^E{R,,iz){\RN{z)\%,) (3.31) 

2 

+ {R,,iz)i\RN{z)\%j) + rN, (3.32) 

where as before rjv contains the third cumulant term corresponding to p = 2 in 
(|7.ip for /c 7^ i, and the error terms due to the truncation of the decoupling formula 
(|7.ip at j3 = 2 for fc 7^ I and at p = for k = i. One can treat as before and 
obtain 

rN { — I , (3.33) 

uniformly on C \ R. Indeed, the bound on the truncation error term at p = 2 for 
k ^ i follows from the fact that it can be written as a sum of 0{N) terms, where 
each term is bounded by O \ 3m zj"^) . The truncation error term at p = 

for fc = i contains one term bounded by O ^ ^''j^^-' | 3mz|~^^ . The third cumulant 
term can be written as 



K3 
2I7V3/2 



E 



Rij{z) 



^ dRkjjz) dRijjz) 
dXik dXik 



Rkj (z) 



(3.34) 
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We will treat each of the three subsums in p.34|) separately. The first one is equal 
to 



1^3 n^rV^ d'^Rkjiz) 



ME 



K3 



- 2|^3/2-[4E($^ i?y(z)i?,fe(z)i?fefe(z)i?y(z)) + 2E(^ R,,{z)Rkk{z)Rkj{z)R,j{z)) 
2-E{Y,iR^k{z)fRok{z)R,,{z))]- 



The same arguments as after (|3.13p bound it by O ^ jv| 3mz|^ ) • second subsum 
equals 



K3 



E 



E 



k:k^i 



dRkjjz) dRijjz) 



K-3 



E 



J2 {Rzj{z)Rkkiz) + Rik{z)Rkj{z)) {Ru{z)Rkj{z) + Rji{z)Rki{z)) 

k:k^i 



It follows from p.l4p that the second subsum is 0( jv| j^zj-^ )■ Finally, the third 
subsum equals 



K3 



2I7V3/2 



E 



E ^kjiz) 

k:k^i 



d^R,,{z) 
dXl 



< 



K3 



2!iV3/2 



Oi\3mz\-')Ej2\Rkjiz)l 



so it is also O y N\3mz\\ 

Using the bound (|3.12p on the variance of tTNRN{z) and p.HI3.2p . we estimate 
the last term in p.30p as 



E{R,,{z)tl■NRNiz)R^J{z)) = gN{z)E{R,,{z)R,,{z)) + O 



1 



ImzA^N 



where we recall that gniz) = EtiNRN{z)- Since the two terms in p.3ip and the 
first term in (I3.32p are bounded by O ^ | jm^iajv ^ , we conclude that 



zE {R,j{z)R,,{z)) = Ei?„ (z)% + CT^gAr(z)E {R,j{z)R,j{z)) + O 
uniformly on C \ K. We now rewrite p.l6p and p.23p as 



n(|3mz|-i) 



N 



zER.,j{z) = S^j + (7^gN{z)ER,j{z) + O 



P4(|3mz|~^) 
N 



(3.35) 



(3.36) 



Multiplying both sides of (|3.36l) by ERij{z) and subtracting from (|3.35p . we obtain 



.2„ ,.^^T..o ..^^_^/'^5(|3"^^r') 



iz-a'gNiz))ViR,,iz))^0 



N 



(3.37) 



uniformly on C \ M. Repeating the arguments after (|3.23p . we conclude that 

'Pe{\3mz\-^y 



YR^J{z) = 



N 
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(3.38) 



This is exactly (|3?3| . 
Proof of (i^ 

Now, we turn our attention to the proof of p.4p . Let us assume that /i has finite 
fifth moment, /ii has finite third moment, and i 7^ j. Without loss of generality, we 
can assume i = 1 and j = 2. We write the Master equation for ERi2{z) as 

zEi?i2(z) =^(j2 ^lkRk2{z)^ = fT^E {Rl2{z)tTNRN{z)) 
2 

+ ^E((i?Ar(2)2)i2) +rjv, (3.39) 

where we apply the decoupling formula (j7.1l) to the term K(XikR2j (z)) and truncate 
it at p = 3 for fc 7^ 1 and at p = 1 for k = I. Thus, the term contains the third 
and fourth cumulant terms (corresponding to p = 2 and p = 3) for fc 7^ 1 as well as 
the error terms due to the truncation of the decoupling formula at p = 3 for fc 7^ 1 
and at p = 1 for fc = 1. We note that in order to truncate the decoupling formula 
at p = 3, we have to require that fi has finite fifth moment. 

It follows from (|3.3p and p.l2p that we can replace in p.38p E {Ri2{z)trNRi\[{z)) 

by E(i?i2(2:))EtrAri?jv(z) up to the error O ( ^'"^?/2' ■ We bound the absolute 

values of the terms in (|3.39p by O — (see Lemmas 13.21 and 1531 below) . 

Combining these results, we obtain 

zERMz) = <7^gNiz)ER,2{z) + O . (3.40) 



Repeating the arguments after (|3.24p . we obtain p.4p . 

Lemma 3.2. Let fi have finite fifth moment and fii has finite third moment. Then 
E((i?^.(z)^)i2)=0(^^|^), (3.41) 

uniformly on C\R. 

Proof. We write the Master equation and use p.2p to obtain 

zE{R\2 ^zEY, RikRk2 = ER12 + E 511] XuRikRk2 

k I k 

= -J^^iy^^jRllRlk + RllRlk)Rk2] + ^E[^ Rlk{RklRl2 + RklRu)] 
l.k l,k 

+ o(^iM)),,„, 

where r^r contains the third and fourth cumulant terms and the error terms due to 
the truncation in the decoupling formula at p = 3 for I 7^ 1 and at p = 1 for 1 = 1. 
Thus, we have 

zE{R'),2 - <j'E[{R%2trNR] + 2^m')i2 + ^E[i?i2 Yl + ^ (M2ll^ 

l.k ^ 

= a'9^{z)E{R'),2 + ^mi2 E Rfk] + o fM2p:!l\ + 

Lk ^ ^ 



■ rN 
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where we used (13. 9p . We note that 



and using p.3p . we have 

^iv«„,./,v,i: . o [M^) o{^).o 

Thus, we arrive at 

(z - a V(^))IE(i?')i2 - O ( ^'^'^T/f'^ ) + (3-42) 

Rather long but straightforward calculations bound in (|3.42[) by O ^ '^°^^^T/2^ — ^ 
We leave the details to the reader. Therefore, we have 

(z - a'g^{z))m'),, = O ^M20:!l^ . (3.43) 

Now consider, as before, On — {z : |2fmz| > LN^^^^}, where the constant L is 
chosen sufficiently large. It follows from p.25p that 

,2^ 1 ^/P7(|3m2|-i)\ ^/P8(|3mz|-i) 



for z G Oat. On the other hand, if | 3m z| < LN^^^'^, then jvi/2|"'jnt^|2 ^ a-nd 



for z ^ Oat. Lemma 13721 is proven. □ 
Lemma 3.3. The term on the r.h.s. of i3. 39\) satisfies 

Proof. First, wc look at the third cumulant terms 



^3 



-[4E(^ RuRikRkk) + 2E(^ RiiRkkRk2) + 2E(^(i?u)2i?2fe)]. (3.46) 



2!iV3/2 ' 

k^l k^l k^l 

To estimate the first subsum in (I3.46p . we write 

|E[i?i2(^i?ifci?fefc)] -E[i?i2]E[^i?ifei?fefc]| < (V(i?i2))'/' I YiY^RikRkk) 

k^l k^l \ k^l 

Taking into account (|3.2p and (|3.14p . we have 
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Therefore, we can bound the first subsum in p.46p by O (^ ^'^^^^/J — • To bound 
the second subsum in (I3.46p . we note that 

^(Ei?iii?fefe)Ei?fe2 = OiP7{\3mz\-^)) (3.47) 

k^l 

by ([321) and ([3T4l) . To bound 

^(V(i?ni?/cfc))'/'(Vi?fc2)i/2, (3.48) 
Ml 

we use (|3.3p and 

V(i?iii?fcfc) < E (|[(i?ii - gN{z)) + gNizMRkk - 9n{z)) + gN{z)] - .g^(^)|') 
= EK^ii - 9N{z)){Rkk - 9n(z)) + gN{z){Rkk - 9n{z)) + gAr(^)(i?ll - 5Ar(z))p. 
Using (|3.3p , we see that 

V(i?ni?..) = O . (3.49) 

The bounds p.49p . p.3p . and p.47p then show that the second subsum in p.46p is 

^ / P7(| Jm^ri) \ 

Finahy, we bound the third subsum in p.46p by using the estimate 
I ^(-Ru)'i?2fc| < |3mzr3. 

Now, let us look at the fourth cumulant terms: 

^[i8E(^ RiiRikRkkRk2) + 6E(^ Rii{RkkfRi2) (3.50) 

k^l k^l 

+ 18E{Y^iRik)^RkkRi2) + 6EiJ2{RikfRk2)]. (3.51) 

fc^l k 

Clearly, 

\J2RiiRikRkkRk2\ <\3mz\-\ 
|^(i?ife)2i?fefei?i2| < \ 3mz\-\ and 
\J2iRikfRk2\ <\Jmz\-\ 

k^l 

To estimate the term -^^{J2k^i Rii{Rkk)^ R12) , we note that by using p.2j|3.3p 
^ E(i?n (i?fefe ) ' )Ei?i2 = O (Ps ( I 3m z r 1 ) ) . 

We are left with estimating 

Y,(y{Rii{Rkkf)Y'^ (Vi?i2)'/' <iV|amzr3 (Vi?i2)'/' = 0(^i/'P6(|amz|-i). 

Combining the estimates of all fourth cumulant terms, we observe that the sums in 
(13.50113. 511) are bounded by O ( 
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To bound the error term, we note that it contains 0{N) terms, such that 
each of them is at most O | j^^pjv5/2 ^ ■ Thus, the error term is bounded by 

^ { \dmz}-m/^ ) ■ Lemma 231 is proven. □ 

This finishes the proof of Proposition 13.11 □ 

4. Expectation and Variance of Matrix Entries Of Regular Functions 

of Wigner Matrices 

In this section, we estimate the mathematical expectation and the variance of 
matrix entries f{XM)ij for regular test functions /. As before, without loss of 
generahty, we can restrict our attention to the real symmetric case. The proofs in 
the Hermitian case are very similar. 

To extend the results of Proposition l3.1l to the case of more general test functions, 
we will exploit the Helffer- Sjostrand functional calculus (see [21], [15]) that depends 
on the use of almost analytic extensions of functions due to Hormander [52] , [53] . 

We start by proving Proposition [5TTJ 

Proof. Let us first assume that / has compact support and prove (|2.19p . Using the 
Helffer-Sjostrand functional calculus (see [H], [H]), we can write for any self-adjoint 
operator X 



(4.1) 



^^''^ = -n id^zV^'^'^'^y ' d'z'-^2[d^+%) 
where: 

i) z = X + iy with a;, y € M; 

ii) /(z) is the extension of the function / defined as follows 

/W:=(i:»^ («) 

Tl = 

here a G C°°(R) is a nonnegative function equal to 1 for \y\ < 1/2 and 
equal to zero for \y\ > 1. 

It should be noted that the r.h.s. of (|4.ip does not depend on the choice of I and 
the cut-off function a{y) in ^i^{see e.g. Theorem 2 in [15]). For X = Xn, (|4T]) 
implies 

f{XNh = -- I %Ru{z)dxdy (4.3) 

TT Oz 

It follows from ([5T]) that with I = 6 in (g^J) 

EfiXNh = -E- / %Ru{z)dxdy (4.4) 

TT Jc OZ 



= [%9Az)dxdy~- [ ^ eii{z)dxdy 



(4.5) 



f{x)d^sc{x) -- [ ^eii{z)dxdy (4.6) 

TT ./r OZ 



where 

|e,,(z)| = \RR,^{z) - g,{z)\ < Ci 



N |/m(z)|6 
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and Ci depends on /i, /ii and L, where supp{f) C [— L, L]. 
Using the definition of / (see (|4.2I) ') one can easily calculate 

df Ifdf .df\ 

— = -I — (4.7) 



i9z 2\dx dy 



(4.8) 



n=0 

and derive the crucial bound 

d~f 



dz 



< 



consi|l/|lc.v([_L,+L])|y|', (4.9) 



for / G C'J(M) with supp{f) C [— Therefore, the second term on the r.h.s. of 
~5l) can be estimated as follows 

I- / %eu{z)dxdy\ (4.10) 

< - [ \% eu{z)\dxdy (4.11) 

< Ciconst{L,fi,fii)\\f\\c7{[-L,+L])-^ J dxxf{x) J dyxaiy) (4.12) 

where Xf ^.nd Xo- are the characteristic functions of the support of / and of a 
respectively. This proves ([^1^ for / e Cj(M). 

To prove (|2.20p . one has to generalize the estimate p.ip in Proposition 13.11 to 
the whole complex plane. We claim the following bound 

\ERu{z) ~ ga{z)\ = \EXx^R^ - 5,(z)| < (|z| + M) ^"^'^]^^' ^\ (4.13) 

uniformly in z G C, 3m z 7^ 0, where M is some positive constant. The bound (I4.13P 
follows from p. 161) . The proof is identical to the proof of similar bounds given in 
[T^ Proposition 4.2 and [TT] Section 3.4. Using the Helffer Sjostrand functional 
calculus as before, one proves (|2.20l) provided / has eight continuous derivatives 

and (|x| + l)^(a;), < Z < 8, are integrable on M. 

In the case of the off-diagonal entries, one takes ^ = 5 in (|4.2p . so 

and proceeds in a similar fashion. The only significant difference is that one has to 
replace the estimate (|4.9I) by 

Qf ( S f \ 

-rtix + iy) <consimax |--4(a;)|, 1 < / < 6 (4.15) 

oz \ dx'- ) 

which together with (|3^ implies ((2?2T|) . To prove (j2.23l) . one takes / = 9 in (14. 2p 
and uses (|3.4p . 

To bound the variance, we write f(XM)ij using (|4.3p with 

n=0 
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Then 



^ = (4.17) 



dz 2\dx dy 

x){iyT Y 

"dy ' 2' ^■"^^"'^^ 3! 



n=0 



and, in particular, 



(x + iy) <Constm&x(\^{x)\, l<l<4:]\y\^ . (4.19) 



9/ f d'- f 

— (x + iy) < Const max ( I ^— J 

Now we are ready to bound Y{f{X]\f)ij). We write 

V(/(XAr)y ) = V l^-i ^ ||i?,^.(z)da;dy^ = i_ £ £ ||||coz;(i?„ (z), i?,, («;))dxd2/dud^;, 

(4.20) 

where z = a; + iy, w = u + iv. We then obtain the following upper bound from (13. 3p 
V(/(X^)y ) <^ J^J^ I li 1 1 1^ I (^)) Vv(i?,, {w))dxdydudv (4.21) 

Const / /■ ,9/, 1 , , , \ ^ 



N \Jc 'dz'\3mz\3' 

Plugging (|4J6|) in (|4:22|) and using (|4.18p . we obtain (|2.22p . Proposition 12.11 is 
proven. □ 

If /i satisfies the Poincare inequality (|2.15|) . one can generalize the results of 
Proposition O Recall that we defined |/|£,r and |/|£,5 in (|2.27ll2.28p . 

Proposition 4.1. Let Xn = "^^^^v be a random real symmetric (Hermitian) 
Wigner matrix with the marginal distributions /i and /ii of the matrix entries sat- 
isfying the Poincare inequality 112. 1 5\) and / : K — > R 6e a Lipschitz continuous 
function on [—2a — S,2a + S] for some S > 0. Let us assume f satisfy the subexpo- 
nential growth condition ^2.37\ l. Then 

P ilfiXNh - E{f{XN),j)\ > t) (4.23) 
< 2Kcxp -irjr— ]+{2K + o(l)) exp — , 



2|/l 



C,5 



where \f\c,5 is defined in h2. 28\) . K is defined in !i2.39\) . and v is the constant in 
the Poincare inequality 112.15]) . If in addition, f G C^(M) for i = j ( f & C^(M) 
for i ^ j), then 

nf{XN\,) - 5,j f{x)^^Aa^-xHx + o(^^y (4.24) 

// the marginal distributions fi and fii satisfy the Poincare inequality \2. 1 5]) and 
f is a Lipschitz continuous function on R, then 

P {\f{XNh - ni{XNh)\ > t) (4.25) 
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where |/|£,r is defined in ^2.27^ . 

Proof. Let us assume that fi and /ii satisfy the Poincare inequahty (|2.15p . Suppose 
that f{x) is a Lipschitz continuous function on M with the Lipschitz constant |/|£,r. 
Then the matrix- valued function f{X) on the space oi N x N real symmetric (Her- 
mitian) matrices is also Lipschitz continuous with respect to the Hilbert-Schmidt 
norm ([T4j. Proposition 4.6, c)) with the same Lipschitz constant. Namely, 

WfiX) - fiY)\\„s < \f\cM\\X - Y\\hs, (4.26) 

where 

\\X-Y\\HS^{M\X-Y\')y^\ (4.27) 

(We note that even though (|4.26p was proven in [T3] only for real symmetric matri- 
ces, the proof for Hermitian matrices is essentially the same). Therefore, f{X)ij is 
a Lipschitz continuous function of the matrix entries of X with the same Lipschitz 
constant. Since the Poincare inequality tensorizes (|19j. [1]), the joint distribution 
of the matrix entries {Xa^ ^ < i < X, \/2Xjk, 1 < j < fc < N} of Xn satisfies the 
Poincare inequality with the constant ■^Nv. Therefore, for any real-valued Lipschitz 
continuous function of the matrix entries with the Lipschitz constant 

\GiX)-G{Y)\ 
---^-^^ WX-YUs ' 
one has (see e.g. [I], Lemma 4.4.3 and Exercise 4.4.5) 




P(|G(Xjv) - EG{Xn)\ >t)< 2if exp | -i^T^t ) , (4.28) 
with 

K = -^2nog(l - 2-^-'). 

i>0 

Applying to f{XN)ij one obtains (IT^ . 

Now let us relax our assumptions on / and consider / : M — > M Lipschitz contin- 
uous function on [— 2(t — (5, 2(t + 5] for some 5 > that satisfies the subexponential 
growth condition (|2.37p . Let h{x) be a C°°(M) function with compact support that 
is identically one in the neighborhood of the support of the Wigner semicircle law, 
i.e., 

h{x) = l for a; e [-2cr-(5,2(T + (5], /i e C;?°(]R). (4.29) 
For non-constant /, we can always choose h in such a way that 

\hf\cM^\f\c^ (4.30) 

Note that 

f{XNh = Uh){XN),, when \\Xn\\ <2a + 5. (4.31) 
It follows from the universality results on the distribution of the largest eigenvalues 
of Xjv (see [24] and also [18], [40], [31], [32], [39], [41], [42]) that 

11X^11 =2fT + 0(7V-i/2-i/i"") 

with probability going to 1. Moreover, G{X) = \\X\\ is a Lipschitz continuous 
function of the matrix entries with Lipschitz constant one. Thus, (|4.28p implies 



\Xn\\ - 2cr| >t)< {2K + o(l)) exp --^—t . (4.32) 
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In particular, 



IXatII >2a + 5)< {2K + o(l)) exp I \ . (4.33) 



Then the estimate (j4.23p for f{XM)ij foUows from the estimate (|4.25p for {hf){XN)ij, 
(|OT]) . and (11331). Finahy, the estimate ^jT/M^ fohows from (l^jg)) and for 
{hf){XN)^J, ^M, and gSl. □ 



Remark 4.1. Let f'^^\x) ^ where z ^ [-2a- d,2a + S] for some S > 0. Th 



en 



If z one has 



dist{z,[~2a-S,2a + S])^- ^^'^^^ 



l/^'^kR^n;^^- (4-35) 



In a similar fashion, for f^'^{x) = — , z,w ^ [—2(7 — S,2a + S]. one has 



\z — w\ 

list{z, [-2a -S, 2a 

For z, w ^ M, one has 



< 2 ^ ^ ^. (4.36) 

mm{dist{z,[-2a-S,2a + 6]),dist{w,[-2a-S,2a + S])f 



|/^^-"^k«<2 . J' ,7' ^. (4.37) 

mm(| Jmz|, I Jmw;!)'^ 

Remark 4.2. Applying the Poincare inequality to Rij(z) one can replace the esti- 
mate iTOj) by 

5. Fluctuations of the Resolvent Entries 



In this section, we prove Theorems 12. 1[ [2721 12. 51 and 12. 61 We start with the proof 
of Theorem O 



Proof of Theorem\MJi As in Section [H we denote by X^"'\ and the 

m X m upper-left corner submatrix of matrices X^^Wn, and R^, where m is a 
fixed positive integer. We denote by x'^"™^ the {N — to) x {N — m) lower-right 
corner submatrix of Xjy, and by 

i?(z)= (z/Ar_„-X(^-"))"\ 

the resolvent of X^^^"^\ We will often drop the dependence on z in the nota- 
tion of i? = if it does not lead to ambiguity. In addition, let us denote by 
x^^\ . . . , a;*^™-' S E^~™ the vectors such that the components of 1 < i < to, 
are given by the last N — m entries of the z-th column of the matrix Xn ■ Finally, 
we will denote by B the (A'^ — to) x m submatrix of Xj\j formed by the vectors 
(columns) a;(*\ 1 < i < ni, and by B* its adjoint matrix. 

Since the fourth moment of /x and the second moment of /xi are finite, ||AAr|| 
converges to 2a almost surely ([6]). Thus, for fixed z G C \ [—2a, 2a], R — (z x 
iN-m — A'-^^™'')^^ exists with probability 1 for all but finitely many N (obviously. 
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R always exists for Jm z ^ 0). Moreover, the mxm upper- left corner of the resolvent 
matrix Rn{z) — {zIn — Xn)^^ , denoted by R^"^^{z), can be written as 



i?('"'(z) = - - = i^zl^ ~ 7^^^"^ " B*RBj . (5.1) 

Let us denote 



T zl^ - - B*RB, so i?'™) 

ViV 



(5.2) 



T = {z- a'^ga{z)) I„, j^^n = —^Im ]=Tn{z), (5.3) 



Write 



where 

{TnMz) = Trj{z) = W^j +Vn (^{x^'\Rx^^'>) - a^g4z)S.,j') , 1 < i,j < m. (5.4) 
We rewrite (|5.4I) as 

r^(z) = vF('") + rjv(^), (5.5) 

where the entries of the matrix Yn{z) are given by 

iYNiz)h, = Y,j{z) = ViV ((a;«,i?(z)x(^')) - <J^gAz)S^j) , 1 < < m. (5.6) 

Remark 5.1. The Central Limit Theorem for random sesquilinear forms (see be- 
low) implies that the entries ofYiq^z), and thus the entries ofT]\[{z) as well, are 
bounded in probability. Recall that a sequence {£,n}n>i of M.'^'^ -dimensional ran- 
dom vectors is bounded in probability if for any e > there exists L{e) that does 
not depend on N such that P(|Cjv| > L{e)) < e for all N > 1. 
Then, 

ViV - ff.(z)/™) = gl{z)T^iz) + O (-^) , (5.7) 



in probability (i.e. the error term multiplied by \/ N is bounded in probability) 



Taking into account ()5.4p , ()5.7p , (|5.5p , and ()5.6p , we can prove the weak conver- 
gence of the finite-dimensional distributions oiT]\f{z), z € C \ [—2a, 2a], defined in 
(|2.ip . to the finite-dimensional distributions of T(z), defined in (|2.6p by proving 
the weak convergence of the finite-dimensional distributions of Yj^{z) to those of 
Y{z), defined by (12.7112.121) . 

To this end, we use the Central Limit Theorem for random sesquilinear forms 
due to Bai and Yao [5^ in the form given by Benaych-Georges, Guionnet, and Maida 
in Theorem 6.4 in |S]. For the convenience of the reader, we give the formulation 
of this theorem in the Appendix (Theorem [TTJ ■ 

Let p be a fixed positive integer and zi, . . . , Zp e C\ [~2a, 2a]. To study the joint 
distribution of the entries ((i?Ar)(0;))i, j, , 1 < I < p, I < ii < ji < rn, it is enough 
to study the distribution of their linear combination. Let a^g\, 1 < s < i < 

TO, 1 < i < p, be arbitrary real numbers and consider 

p 

TW^^'*' = (^a'^l mc{R{z,)) + b'^l 3m(i?(z,))) , 1 < s < t < to, (5.8) 
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where for any linear operator A 

<H.(A) = ii±iil, 

Now, we show that the results of Propositions 13.11 and 12.11 and the almost sure 
convergence of ||^Af|| to 2a imply that the conditions (|7.11l [7!T^ of Theorem 17.11 
are satisfied. First, we note that as N ^ oo, 

trjv (me{R{z)) y\e{R{w))J V++{z, w), (5.9) 

ti-AT (^3m(^(z))3m(^(w))^ ip (5.10) 

iiN {^t{R{z)'3m(R(w)^ -^^+_(z,u)), (5.11) 
1 ^ 

- Y,{^z{R{z)))u{^z{R{w)))u ^ ^t{9.{z)) ^t{9.{w)), (5.12) 



1=1 

N 



^ Y^{'3m{R{z))),,{3m{R{z)))u -> 3m(5<,(z)) 3m{ga{w)), (5.13) 

1 ^ 

- Y,{mz{R{z)))u{3m{R{w)))u ^ 9ie(5,(z)) 3m{g,{w)), (5.14) 



i=l 

for z, w G C \ [—2(7, 2cr], where iy9++(z, w), Lp (z, w), and (p^ (z, li;) are defined in 

(|2.3ll2.5p . and the convergence is in probability. 

Indeed, (|5. 9115. 11]) follow from the semicircle law (and the convergence can be 
taken to be almost sure). In particular, for real z and ui, in order to avoid singular- 
ities, one can replace R{z), R{w) by h{X]\f)R{z), h{XN)R{w), where h is defined 
in (|4.29p . and use the fact that II-^^atH — 2(t almost surely as — cx). 

Let us now prove (|5.12p . The proofs of (|5. 13115. 14)) are similar. We can assume 
that 3m z ^ 0, 3mw ^ 0. Otherwise, one has to replace R{z) by h{XN)R{z). We 
write 

\imciR{z)))u{mR{w))h - ^e(g.(z)) V\c{g4w))\ 

< |(5He(i?(z)))., - mcig4z))\\{mciRiw)))u\ + \ mc{gAz))mc{R{w))h ~ ^e(g.(u;))| 

< \{D\z{R{z)))u - mcig.{z))\—^ + -}—\{D\c{R{w)))u - ^t{ga{w))\. 



Thus, it follows from ()3.3p that 

E|($Ke(i?(z))),,(fHe(i?(u;))),,-9ie(5.(z))^ne(g,(«;))| < ( + ^r^) 0{N-'/'), 

\ I Jmwl I Jmz\J 

which implies (|5.12p . 

It should be noted that Theorem l7.1l is proven for non-random matrices A^^'*"* , 1 < 
s < t < m. Since the convergence in probability does not imply almost sure con- 
vergence, an additional argument is in order. Let A^^'*^ be defined as in (|5.8p . and 
— I < i < m. By calculating the second moments of 
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one can show that these random variables are bounded in probabihty. Therefore, 
it is enough to prove convergence for a subsequence Nn — >■ oo. Since convergence in 
probabihty imphes almost sure convergence for a subsequence, we can now apply 
Theorem 1 7 . 1 1 directlv to a subsequence. 

Applying Theorem I7.1[ we establish the convergence of the finite-dimensional 
distributions of yAr(z) and obtain (|2.7II2.12|) . Theorem l2.1l is proven. □ 

The proof of Thcorem l2.5l is very similar to the proof of Theorem 12. 1[ Theorem 
17.11 plavs the central role in our arguments again. We choose 

p 

M%''^ = (a^l y\z{R{zi)) + h^l 3m(i?(zi))) , 1 < s < i < m, (5.15) 

i=l 

where a\ i, b), [, 1 < s < t < m, 1 < i < p, are arbitrary complex numbers and 

ai*s, 6i*s, ^ ^ s < m, 1 < i < p are arbitrary real numbers. Applying Theorem 
6.1, we establish the convergence of the finite-dimensional distributions of Yn{z) 
and obtain (I2.42ll2.47p . 

Now, we prove Theorem 12.21 

Proof of Theorem \2.2i Let /i and fii satisfy the Poincare inequality (I2.15p . To prove 
the functional limit theorem, i.e. the convergence in distribution of the sequence 
of probability measures Vn on C(X',C'"(™+i)/2), it is now sufficient to prove that 
the sequence Vn is tight ([10]). For this, we need to show that for 

Tjv(z) - Vn (i?(™)(z) - g,(z)/™) , zeC\ h2cT, 2a]. 

the following conditions are satisfied: 

(a) for some fixed zq e V, for each e > 0, there exist sufficiently large K and 
A^o such that 

P(||T7v(zo)|| > if) < e foranA^>Afo, (5.16) 

and 

(b) for each e > and a > there exist 7 > and Nq such that 

P(wT„ (7) > a) < e for aU N > Nq, (5.17) 
where a;Tjv(7) denotes the modulus of continuity of Yn{z) on V, namely 

'^Tjv(7) = swp|^„^|<^||Tjv(z) - Tjv(w)||, (5.18) 

where the supremum is taken over all z,w ^ V such that — if | < 7. 

The property (a) immediately follows from the definition of Tjv and the bounds 
p.m3.3p in Proposition l3.1l To prove (b), we replace Rn{z) by h[XN)RN{z) in the 
definition of T n{z) (if I? n M = 0, this procedure is not needed), where h is defined 
in (|4.29l) in such a way that supp{h) n I? = 0. We note that Rn{z) = h{XN)RN{z) 
almost surely for all z and for sufficiently large N. 

It then follows from the results of Proposition 14. II that uniformly in z,w V, 

P(|TAr(z) - Tn(w)\ >t)< 2Ke^p ( -const i- — - — - ) +(2i4:+o(l)) exp 

V \z-w\J 

(5.19) 
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with some consti > 0. In addition, for any z d V, 

di^Nhiz)^ > A < 2Kexpi-const2t) + i2K + o{l))exp ( ) ^ (5.20) 



dz 

for some const2 > 0. Without loss of generahty, we can assume that I? is a rectangle 
with the sides parallel to the coordinate axes. We then partition T) into 0(2^") 
small squares UiDj^"-* of the diameter 2~", l<n < ni — constlog{N){l + o{l)), 
where const > is chosen so that 

const > (5.21) 
2 ^ ^ 

We then estimate the probability of the event that 

s„ sup ||Tjv(2) - TAr(w)|| > Ani™2-", (5.22) 

where the supremum in (j5.22l) is taken over all pairs {z,w) that are the vertices of 
the same small square. Using (|5.19p . one can show that this probability is 



constiA iQQ,\ , ^ / / \/vN5 ^ 



O (^exp( "-^n'"")j + O \^cxp{- 

uniformly in N. We can also estimate the probability of the event 

:^sup|l^^^^J^|| > Anr-^constioO(bg(iV))ioO(l + o(l)), (5.23) 
where the supremum in (|5.23p is taken over all vertices of the partition UiD^^^^\ 

by 

/ const2A loo^^ , n I f y^^NS \ 
O exp( n^'") +0 exp( — ) , 



uniformly in N. 

Finally, we note that since we choose const > in ni = constlog{N){l + o(l)) 
to be sufficiently large so that (|5.2ip is satisfied, we have 

|^„^_sup||^(I^|||<l, (5.24) 

since the second derivatives of the entries of Tn{z) are trivially bounded by 
const^y/N , where consts depends on 7). Now choosing A sufficiently large, we can 
make the probability P (cjy„ (ca| log a\) > a) smaller than e for a suitable constant 
c > 0. We leave the details to the reader. □ 

6. Fluctuations of Matrix Entries of Regular Functions of Wigner 

Matrices 

We give the proofs in the real symmetric case (Theorems 12.31 and 12.41) . The 
proofs in the Hermitian case (Theorems 12.71 and 12. 8p are very similar. First, we 
assume that /i and /ii satisfy the Poincare inequality (|2.15p and prove Theorem l2.4l 
Then we will extend it to the case of finite fourth moment and prove Theorem 12. 31 

Proof of Theorem \2.4\ We start by considering a test function / which is analytic 
in a neighborhood of [—2a, 2a] and takes real values on R. We write 

/(^jv).j = ^ / f{z){RN{z)hdz, (6.1) 
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where 7 is a clockwise-oriented contour in the domain of analyticity of / that 
encircles the interval [— 2cr, 2cr]. Then 



(^fiXNh - J ^ f{x)^Vla^^dx^ "2^/ fi-^)0:N{z)hdz, l<i<j<rr 

(6.2) 

By the functional convergence in Theorem l2.21 the r.h.s. in (j6.2p converges in distri- 
bution to the distribution of independent (up from the diagonal) random variables 

^ I /(z)T.,(z)dz= I f(z)gl(z)dz^W,,+^^ I f{z)gl{z)Y,,{z)dz, 

(6.3) 

^ < i < j < Tn. We evaluate 

-L f fiz)gUz)dz ^± f fizf-^^^^dz = (6.4) 

^2^ y^/(-)-5.(^)d- = ^y_^^-/(^)2^v/4a^dx. 

The last integral in ()6.3p is a real Gaussian random vector with independent entries 
and variances 

(6.5) 

Let us first consider the off-diagonal case i 7^ j. By (|2.10II2.1^ and (12. 2112. 5p . we 

have 

E {gl{z)gl{w)Y,,{z)Y,,(ui)) - E {gl{z)Y,,{z)) E {gl{w)Y,,{w)) = a^gl{z)gl{w)^{z,w) 

4 2/ N 2/ .9<y{z) - ga{w) i ( ^9l{z)ga{w) - ga{z)gl{w) 

= 9Az)9Aw) = -cr g„{z)g^{w)- 



z — w z — w 

-a'^g„{z)g„{w) 



2 ^„^„ ,^„^zg^{z)ga{w) - g„{w) - wg„{z)ga{w) + g^{z) 



z — w 

= -a- g„[z)g„{w) ga{z)g„{w) H 

\ z — w 

= -(^'^9liz)9l{w) + {if{z,w) - ga{z)g„{w)) . 

We note that ^p{z,w) — g„{z)g„{w) — Cov{^^, ^^j.^ ), where rj is distributed 
according to the semicircle law (|1.7I) . Then 

J ^ J f(^)f('^)('^(^^'^)~9<T{z)ga{w))dzdw^ 
^ [ ^ [ f{z)fMCov{^,^—)dzdw=Y{f{v)). 

This together with (I0)) . ((0)) proves (P^ . (lOil) for analytic functions. 
In the diagonal case i = j, the previously studied terms contribute to 

E{gl{z)gl{w)Y,,{z)Y,,{w)) - E{g^{z)Y,,{z))E{g,{w)Y,,{w)) 
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with a factor of two. In addition, ()2.7ll2.9p provide one more term 
i^i{fJ-)9Uz)9l{w)- Evaluating 

-L f f[z)gl{z)dz = -L / /(,) (f!^^!)|^(fl^dz = (6.6) 

1 1 o o, , . . 1 



40 ■ / f{z){z'~<j')ga{z)dz^^ r {x^ - a')f{x)-^Vi<J^ - x^dx, 



-2(7 



we prove p.30p . p.33p in the analytic case with the centralizing constant /^jo- fi^) ^ 

duscidx). 

It follows from (|4.32p that if / satisfies the subexponential growth condition 
p.37p on the real line then we can choose the centralizing constants to be E/(Xjv)ij 
in ()2.30p , ()2.33p . To extend the results of Theorem 12.31 to a more general class of 
functions we apply a standard approximation procedure. If / is Lipschitz continuous 
on [—2(7 — (5, 2(7 + S] and satisfies (|2.37p . we choose a sequence of analytic functions 
{fn}, n > 1, such that 

/n(0) = /(O), n > 1, and |/„ - f\c,s ^ as rn> oo, 

where \f\c,s has been defined in p.28p . Let us also choose h in such a way that 
/i : M — )■ M is a smooth function with compact support, h{x) = 1 for a; G [—2 x 
(7 - (5/2, 2(7 + (5/2], and h{x) = for > 2a + jS. We observe that for any n > 1 

/(Xn) + (/M(Xjv), /„(Xjv) + {fnh)iXN), 

with probability exponentially small in \/N. In addition, 

E|(/(l - h))iXNhj\ = 0{expi-constVN)), 

E|(/„(l - h)){XN)^J\ = 0{exp{-constnVN)), n > 1, 

where const > 0, constn > 0. We then choose n sufficiently large so that 
\fh—fnh\c < £■ As in the proof of Proposition l4. 11 we use the fact that for any Lips- 
chitz continuous /, the function f{X)ij is a Lipschitz continuous function of the ma- 
trix entries of X. Therefore, we can show that (^\/N{{fh){XN)ij — {fnh){XN)ij)^ 

can be made arbitrary small (uniformly in N ) for sufficiently large n if we apply 
the concentration inequality (|4.28p to {fh){Xt^)ij — {fnh){Xfq)ij. Finally, we ob- 
serve that uj\fr,) ^ uj^U), a^ifn) ^ a^if), P^{fn) ^ /32(/), d\fr,) ^ d^U) as 
n — )■ oo. 

It follows from Proposition 12.11 and (|4.32p that if / is seven times continuously 
differentiable on [—2(7 — (5, 2(7 + 5] (six times continuously differentiable on [—2 x 
a — 5,2(7 + S\ in the off-diagonal case i ^ j) and satisfies p.37p then one can replace 
E/(XAr)jj by (5y / f{x)dnsc{dx) in (1^311)) . (^351) since 

E/(Xjv)., = (5y j f{x)dfXsc{dx) + 0(1). 
Theorem 12.41 is proven. □ 



Now, we prove Theorem 12.31 assuming only that /i and fii have finite fourth 
moments. The role of (14281) will be played by ([2221) ■ 
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Proof of Theorem \2.3i By Theorem 12 . 1 1 and Proposition 13. ![ we have the result for 
finite hncar combinations 



k 

fix) =\2ai^—, zi ^ [-2(7, 2a], 1 < / < fc, (6.7) 

^ Zl - X 



1=1 

and, more generally, for 

fe 

fix) ^y^aihiix) , Zl ^ [-2(7,2(7], 1 < ; < A:, (6.8) 

where hi e (M), 1 <l < k, satisfy (I4.29p . Applying the Stone- Weierstrass theo- 
rem (see e.g [3S]), one can show that such functions are dense in C^(K). Therefore, 
we can approximate an arbitrary / e C^(R) by such functions /i(a;)/„(x) in such a 
way that 

suppif) c [~A,A], suppihfn) C [-A,A], 

for all n and sufficiently large A > 0, and ||/ — ^/ri||c''([-yi.yi]) ^ 0, as n — > oo. It 
follows from (I2.22p that V(/(Xjv)y — (/i(x)/„)(Xjv)ij) can be made arbitrary small 
uniformly in n provided we choose n to be sufficiently large. Since a;^(/i/„) — 
w'(/), a\hfr.) ^ ^2(/j^^) ^ ^2(^)^ ^2(;^j^) ^ ^2(j) as n ^ OO, the 

result follows. Theorem 12.31 is proven. □ 

7. Appendix 

In our analysis, we need to study the expectation of the random matrix entries 
multiplied by functions of the random matrix. In order to handle this we use the 
following decoupling formula '26': 

Given ^, a real-valued random variable with p + 2 finite moments, and a 
function from C — > M with p + 1 continuous and bounded derivatives then: 

p 



E(?</'(0) = E ^is('^'"ne)) + e (7.1) 

^ — ' a! 



a=0 

Where the cumulants of |e| < Csup^ \4'P+^\t)\&i\£,\P+^) , C depends only 

on p. 

For any two Hermitian matrices Xi and X2 and non-real z we have the resolvent 
identity: 

izl - X^)-^ = izl - Xi)'' - izl - Xi)-\Xi - X2)izl - X2)-' (7.2) 

If X is a real symmetric matrix with resolvent R then the derivative of Rki with 
respect to Xpq, for p ^ q is given by 

dRki 

TTT;: — = RkpRql + RkqRpl- (7.3) 
OApq 



If p = q then the derivative is: 

ORm 



dXpp 



— RkpRpi- (7.4) 
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In a similar way, if X is a Hermitian matrix then the derivative of Rki with 
respect to DitXpq, JmXpq, for p ^ q are given by 

dRki 



d9\cXpg 



— RkpRql + RkqRpl, (7.5) 



dR 

: I [RkpRql ~ RkqRpl) ■ (7.6) 



d 3m Xpq 



When p = q then 



— RkpRpi- (7.7) 



dXpp 



We will use the following bounds on the resolvent: 

where by Sp{X) we denote the spectrum of a real symmetric (Hermitian) matrix 
X. (TTSl) implies 

\\RN{z)\\<\3miz)\-' (7.9) 

which also implies all the entries of the resolvent are bounded by | JTn(z)|~^. Sim- 
ilarly, we have the following bound for the Stieltjes transform, g{z), of any proba- 
bility measure: 

\g{z)\<\3m{z)\-' (7.10) 
Below, we state the Central Limit Theorem for random sesquilinear forms of Bai 
and Yao in the form given in [S]. 

(s t) 

Theorem 7.1. Let us fix m > 1 and let, for each N, M}^" , I < s,t < m, be 
a family of N x N real (resp. complex) matrices such that for all s,t, A^^'^^ = 
^A^^'*^^ and such that for all s,t = 1, . . . ,m, 



trN 



^(7,%, asn-^oo, (7.11) 

N 



^El(>'i^''^)-I'^7s, asn^oo. (7.12) 

i=l 

Lei he a sequence of i.i.d. random vectors in^'^ (resp. C.-^ ) such that 

the N coordinates of u^^^ are i.i.d. centered real (resp. complex) centered random 
variables distributed according to a probability measure with variance one and finite 
fourth moment. In the complex case, we also assume that real and imaginary parts 
of each coordinate of u^^-* are independent and identically distributed according to 
a probability measure v on the real line. 

For each N, define the m x m random matrix 

GN:=(^VN(^{u^P\M^^''^u^^^)~6p,tr^{M^^^P^))), l<p,q<m. (7.13) 

Then the distribution ofG^ converges weakly to the distribution of a real symmetric 
(resp. Hermitian) random matrix G = (gp.q), 1 l£ P, q l£ m, such that the random 
variables 

{gp,q, ^ < P,q < m}, (resp.{gs,s, I < s < m, ^cgp^q, 3mgp^q, 1 < p,q < m}) 
are independent for all s, g^s ~ ^(0, 2a^ ^ + ^4(2^)7^) (resp. 
gss ~ N{0, al ,. + iK4(t/)7s)), and for all p ^ q, gs,t ^ N{0, a^ g), 
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(resp.) d\tgs,t ^ iV(0, ^(Tp 3mgs,t ^ -^(0, ^cr^.g), where Ki{v) denotes the 
fourth cumulant of v. 

Remark 7.1. Almost simultaneously with our paper, L.Pastur and A. Lytova 
posted a preprint 134] where they extended the technique of [28 and gave another 
proof the convergence in distribution for a normalized diagonal entry \/iV'(/(Xjv)ii — 
E(/(XAr)ii) under the conditions that a real symmetric Wigner matrix Xj^ has i.i.d. 
entries up from the diagonal and the cumulant generating function log(Ee^^^^ ) is 
entire. Pastur and Lytova require that a test function f satisfies 

^(l + |fc|)3|/(fc)|dfc<0O, 

where f{k) is the Fourier transform of f. 
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